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Abstract. We consider propagation, storing and retrieval of slow light (probe beam) 
in a resonant atomic medium illuminated by two control laser beams of larger intensity. 
The probe and two control beams act on atoms in a tripod configuration of the 
light-matter coupling. The first control beam is allowed to have an orbital angular 
momentum (0AM). Application of the second vortex- free control laser ensures the 
adiabatic (lossles) propagation of the probe beam at the vortex core where the intensity 
of the first control laser goes to zero. Storing and release of the probe beam is 
accomplished by switching off and on the control laser beams leading to the transfer 
of the optical vortex from the first control beam to the regenerated probe field. A part 
of the stored probe beam remains frozen in the medium in the form of atomic spin 
excitations, the number of which increases with increasing the intensity of the second 
control laser. We analyse such losses in the regenerated probe beam and provide 
conditions for the optical vortex of the control beam to be transferred efficiently to the 
restored probe beam. 
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1. Introduction 

Over the last decade there has been a great deal of activities in slow [H El [3l IH E] , stored 
illTllHlElttOlElinittSllIllIianniinittH] and stationary [HI gOl [HI [221 1231 122 [25] light. 
It was demonstated that a resonant weak pulse of light (to be referred to as the probe 
light) can propagate as slow as several of tens of meters per second [1] in an atomic 
medium driven by a stronger (control) laser beam. The application of the control laser 
beam makes the resonant and opaque medium transparent for the probe beam due the 
electromagnetically induced transparence (EIT) [261 [271 [281 [291 [30] , both beams making 
a A configuration of the atom-light coupling. The EIT can be used not only to slow down 
dramatically the hght pulses [B [21 [SI [2 [5] , but also to store them [3 [SI [131 [ig [161 [IH [T8] 
in atomic gases. The storage and release of the probe pulses has been accomplished 
[3 [El [iHl [IHl [HI [TTl [IE] by switching off and on the control laser [6j . The slow and storred 
light can be applied to reversible quantum memories [H [H [TOl [131 [211 [301 [211 [321 [33] and, 
in the case of moving media [311 [SHI [33 [33 [3E1 [321 [301 [S] > also to rotational sensing 
devices. 

The orbital angular momentum (0AM) [42l [43] provides additional possibilities in 
manipulating the slow light. The optical 0AM represents a new degree of freedom which 
can be exploited in the quantum computation and quantum information storage [13]. 
Most of the previous studies on the vortex slow light have concentrated on situations in 
which the incident probe beam carries an 0AM [HI [l5l [HI [l6l W7\ . 

Here we consider another situation in which the incident probe beam does not 
have an optical vortex, yet it gains the 0AM when retrieved. For this purpose one 
of the control laser beam is assumed to have an optical vortex during either the 
storage or retrieval of the slow light. The intensity of such a control beam goes 
to zero at the vortex core, potentially leading to the absorption losses of the probe 
beam in this area. To avoid the losses an extra control laser without an optical 
vortex is used. This makes a more complex tripod scheme of the atom-light coupling 
[l8l[l9l[50l[53[5a[52[Ml[55l[56l[53[58l[5^ The total intensity of the control lasers 
is then non-zero at the vortex core of the first control laser preventing the absorption 
losses. Subsequently we analyze additional losses appearing because a part of the stored 
probe beam remains frozen in the atomic cloud in the form of spin excitations during 
the exchange of the optical vortex between the control laser and the regenerated slow 
light. A number of such frozen excitations increases with increasing the intensity of the 
second control laser. We provide conditions for the optical vortex of the control beam 
to be transferred efficiently to the restored probe beam. 

2. Equations for the probe beam 

In this Section we shall derive general equations for the propagation of the probe beam 
in an ensemble of tripod atoms. In doing so we shall not make use of the usual slowly- 
varying amplitude approximation [55l [53 [56] for the spatial coordinates in the direction 
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Figure 1. Light induced transitions in the cloud of cold atoms involving a tripod 
configuration of the internal atomic states. 



perpendicular to the light propagation. This will make it possible to analyse situations 
where the control and/or probe beams have the 0AM. 

2.1. Coupled equations for the probe beam and atoms 

Let us consider an ensemble of tripod-type atoms characterised by three hyperfine 
ground levels 1, 2, and 3, as well as an electronic excited level (figure H]). The 
translational motion of atoms is represented by a four component column operator \I/(t*), 
where the components \l/i(r,t), \E'2(r,t), \l/3(r,t), and 'ifo{r,t) are the field operators 
describing to the center of mass motion in the four internal atomic states. The quantum 
nature of the atoms comprising the medium will determine whether these field operators 
obey Bose-Einstein or Fermi-Dirac commutation relations. The atoms interact with 
three light fields in a tripod configuration of the atom-light coupling [HI [531 ESI [561 ISZ] • 
Specifically, two strong classical control lasers induce transitions |2) — )■ |0), |3) — t- |0) 
and a weaker quantum probe field drives a transition |1) — |0), as shown in figure [1] 
The electric field of the probe beam is 



where u = ck is the central frequency of the probe photons, k = zk is the wave 
vector, and e±z is the unit polarization vector. This field can be considered to be 
either a quantum operator or a classical variable. We have chosen the dimensions of the 
electric field amplitude £ to be such that its squared modulus represents the number 
density of probe photons. The probe field is assumed to be quasi-monochromatic, so 
the amplitude S{r, t) changes little over the time of an optical cycle and thus obeys the 
following equation: 



where have introduce the slowly varying atomic field operators $1 = \E'ie^'^^*, $2 = 
^^giK+a;-a;,2)t^ ^ ^^gi(L^i+^-a;,3)t ^ = ^j^j^ f^^^ l^gj^g ^^^e energy of 




(1) 




(2) 
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the atomic ground state 1, ujc2 and oOcs being the frequencies of the control fields. The 
parameter g = ii^ix) jle^h featured in equation ([2]) characterizes the strength of coupling 
of the probe field with the atoms , [i being the dipole moment of the atomic transition 
|1) — )• |0). Note that, unlike in the usual treatment of slow light, we have retained the 
second-order derivative in the equation of motion ([2]). This allows one to account 
for the fast changes of £^ in a direction perpendicular to the wave vector fc, i.e., in 
the xy plane. Therefore our analysis can be applied to the twisted beams of light 
S{r,t) ~ exp(i/y9) carrying an 0AM hi per photon. 

In the following we shall make use of a semi-classical picture in which both the 
electromagnetic and matter field operators are replaced by c numbers. The equations 



for the matter fields are then 

idt^i = -(7^*$o (3) 

idt% = (woi - i7)*o - fic2$2 - fics'^'a - g^^^i (4) 

idt^2 = W21$2 - ^*c2% (5) 

idt<^-i = c^si'^'a - ^cs'^o , (6) 



Here and fics are the Rabi frequencies of control lasers driving the transitions 
1 2) — )• |0) and |3) |0); 7 is the decay rate of the excited level; 1^21 = 0U2 — uJi + ujc2 — oj 
and = oj-i — (jJi + — oj are frequencies of the electronic detuning from the two- 
photon resonances, wqi = ujq — uji — u \s the frequency of the electronic detuning from 
the one-photon resonance. In writing equations ([S])-© we used the rotating wave 
approximation for the atom-light coupling, and neglected the terms containing atomic 
mass m, since from now on we are not interested in the effects due to the atomic motion: 
The terms containing atomic mass m would be important in equations ([3])-([6]) for the 
description of the light-dragging effects [Ml EZl [39l HO, HI]- Furthermore we assumed 
that the trapping potentials for an atom in the electronic state j {j = 1,2,3,0) are 
position independent and hence can be omitted. 

2.2. Equation for the of probe beam 

To analyze the atomic dynamics it is convenient to introduce the bright and dark states 
of the atomic center of mass motion 

$B = fi-^(l^c2<f2+^^c3$3), (7) 

= n:\n:,<!>2 - , (s) 

where 

= Vl^ufTJn^ (9) 

is the total Rabi frequency. It is the wave-function $b which is featured in the equation 
of motion (jlj) for the excited state wave-function. On the other hand, the dark state 
is not coupled directly to the excited level 0. The original fields $2 and $3 can be 
expressed through <I>b and $d as 

^2 = ^c'i^:2^B + ^c3^D), (10) 
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<I>, = Q^\Q:,^b-^c2<^d). (11) 

Initially the atoms are in the ground level 1 and the Rabi frequency of the probe field 
is considered to be much smaller than Qc- Consequently one can neglect the last term 
in equation ([3]) that causes depletion of the ground level 1, giving idt^i = 0. If the 
atoms in the internal ground-state 1 form a BEC, its wave function $i = ^/nexp{iSl) 
represents an incident variable determined by the atomic density n and the condensate 
phase Si. The latter phase will not play an important role in our subsequent analysis, 
since we are not interested in the effects due to the condensate dynamics. 

Suppose the control and probe beams are tuned close to the two-photon resonance. 
Application of such laser beams cause electromagnetically induced transparency (EIT) 
in which the transitions |2) — )■ |0), |3) — )■ |0), and |1) |0) interfere destructively 
preventing population of the excited state 0. The adiabatic approximation is obtained 
neglecting the excited state population in equation (jlj), giving 

<^B = -9^^S. (12) 

When the adiabatic approximation is valid, the bright and dark states are coupled 
weakly. Therefore, combining equations ([5])), (E]) and neglecting the dark state 
contribution we arrive at the equation for the bright state wave- function 

i5t$B = - ^^c$0 , (13) 

where 

6 = a;2i|6l' + ^31161' - K^2dtC2 + ^3^*^) (14) 
is the two-photon frequency mismatch, with 

^2 = n,2/n, , ^3 = ^cs/^c . (15) 

The derivation of equation (fT3!) can be found in more details in our earlier work [60] on 
the light-induced gauge potentials for the A type atoms. Equation ( fT3ll relates $o to 
the bright state as 

% = n^\-idt + 6)^B. (16) 

Finally, equations ([2]), ( lT2i) and ( |T6l) provide a closed equation for the electric field 
amplitude £: 



This equation applies a wide variety of phenomena. In particular it can be used to 
model light storage by introducing time dependence in or light dragging due to 
spatial variation of Qc or $i. 
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2.3. Non-adiabatic corrections 

In the next Section we shall consider the releasing of the stored light. For this we should 
include non-adiabatic corrections to the equation of motion ( |T7|) . This can be done in 
the following way. From equation (j4]) expressing the bright state $b and substituting 
equation (fT6l) for $o we get 



$B = + ^:\-idt + ujoi - h)n;\-idt + 5)$B . (18) 

The term with the decay rate 7 is larger than other non-adiabatic corrections in the 
above equation. Keeping in the non-adiabatic corrections only the terms proportional 
to the decay rate 7 and neglecting, for simplicity, the two-photon detuning S we obtain 

$B ~ - ^9i$B . (19) 

The solution of this equation, assuming that the control beam is switched on suddenly 
at t = an then changes slowly during the characteristic relaxation time j/^l^ is 

$B = <J>B(0)e-?* - |^<l>i (1 - e-?*") . (20) 



Using equations ( 1161) and ( ]20l) . the equation for the electric field ([2]) takes the form 



+i^(-idt + 5)<^BiO)e-^'. (21) 



At t = the probe field is off, and the information on the previously stored probe beam 
being contained in the atomic coherence $b(0). The regeneration of the probe beam is 
described by the second term on the r.h.s. of equation (12T!) representing the source for 
the electric field. Retaining only the temporal derivatives in equation fl2T]) . we get the 
equation describing the generation of the electric field: 

^ "'■^^9i^ = -^[r?$i^ + (]c$B(0)]e-?*, (22) 

7 

with the initial condition £{0) = at t = 0. For time in access of the relaxation time 
7/fi^ the regenerated probe field evolves to a steady-state value complying with the 
adiabatic condition (|T2l) 

S = -^^b{0). (23) 

In this way, the regenerated electric field S is indeed determined by the initial atomic 
coherence $b(0). The subsequent evolution of the probe field is described by the 
adiabatic equation of motion 0171) containing both the temporal and spatial derivatives 
subject to the initial condition fl23l) . 
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2.4- Co-propagating control and probe beams 

Suppose the probe and control beams co-propag ate: S{r,t) = 8{r,t)e'^', ^c2{r,t) = 
Q2{r,t)e^''''^^ , Qcsi^jt) = Q3{r,t)e^^''^^, where kc2 and fccs are the wave numbers of the 
control beams. For paraxial beams the amplitudes £{r,t), Q2{i^,t) and fi3(r,t) depend 
weakly in the propagation direction z compared to the variation of the exponential 
factors. Equation f|T7j) for the probe field takes then the form 



dtS + v„ 



where we have replaced by its transverse part V^^ = d'^ jdx^ + &^ jdy^ because of the 
paraxial approximation. Here 

2 \ -1 

i+y (25) 

is the radiative group velocity. The term with spatial derivative djdz in equation ( l24l) 
describes the radiative propagation along the z axis with the group velocity fg. 

3. Storing and releasing the light 

The probe beam E'^^^ enters an atomic medium which is illuminated by two control 
beams characterized by the Rabi frequencies V!^^2 ^cs > where the index (s) refer 
to the stage of storing the light. At the boundary the probe beam is converted into 
a polariton propagating slowly in the medium with the velocity vf^ <^ c. At certain 
time t = t^^^ the whole probe pulse enters the atomic medium and is contained in it. 
Since the atomic population is created exclusively by the incident light field, the atomic 
dark-state $d is not populated and, according to equation (fT2|) . the bright-state is 

^S{t^'^) = -9^^'^''\t(^^). (26) 

i Lc 

Equations ffTU]) and ffTT]) give the atomic fields: 

<^>? = , ^? = 6t^i^ ■ (27) 

To store the slow light, both control fields are switched off at t = t^^"^ in such a way that 
the ratios ^^2'' ^c3 remain constant whereas Q^c"^ — > 0. The stored atomic coherences 
no longer have the radiative group velocity and thus are trapped in the medium. To 
restore the slow light propagation, the control fields are switched on again at t = t^^^ 

(v) (r) 

with relative Rabi frequencies ^^2 • The latter can differ from the original ones 

(s) (s) 

Q2 and (^^3 , so the dark-state $d can now be populated. Shortly after the beginnig of 
the release of hght (at t = t^^^) the generated electric field reaches a steady state value, 
as described in the Subsection 12. 3[ Equation fl23|) yields the restored probe field: 

^«(tW) = -^(d?df + &'&^i\t^'^) , (28) 
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where equations (I7l), ( fTSl) and ( l27l) were used to relate the bright state ^^■'(t^''-*) of the 
restoring stage to the stored one ^^''(t^^'*) . 

Since a typical length of the atomic cloud is not much larger than the length of the 
laser pulse, we will assume that the stored and restored electric fields do not change 
significantly during their propagation inside of the atomic cloud. Furthermore we will 
assume that the control beams are abruptly switched off during the storing stage and 
then switched on in the same way when they are restored. Using equations (fT2|) and 
( |28|) . the restored field may be written as: 

^^'' = %)i^cki^*+^ckth^^'^- (29) 

3.1. Control beams with the same spatial behaviour 

Supposed first that the Rabi frequencies of the restored control beams are proportional 

frl is) (r) is) 

to the original ones: Vt).2 = aVt).2 and Vt).^ = af^cs • For slow light this implies that 
d? = d? = d?- Since ^c2C2 + UCs = 1, equation ([29]) yields the following result 
for the regenerated electric field: 

£«(t«) = -^$(^)(t(^)). (30) 

The above relationship represents the initial condition for the subsequent propagation 
of the probe beam £ governed by the equation of motion ( 12^ . The regenerated electric 
field is seen to acquire the phase from the bright polariton at its storage stage and the 
amplitude of 8 is modulated according to at the release stage. 



3.2. Transfer of optical vortex at the retrieval of the probe beam 

(s) 

Suppose now that initially we have a A system with only a single control field: = 
and hence 1^2'' I = 1- Oii the other hand, a tripod system is used in the retrieval stage 
where generally both ^ind are non-zero. In that case equation ( p8|) provides the 
following result for the regenerated electric field: 

£^'\t^'^) = -^&^i\t^'^). (31) 

The equation (1311) represents the initial condition for the subsequent propagation of the 
probe beam S governed by the equation of motion fl2^ in the medium. 



If the second control beam has an optical vortex at the restoring stage f2^2^ ~ e^^*^, 
the regenerated electric field £^^^ ~ e^^"^ acquires the same phase, as one can see from 
equation (PT]) . This means the restoring control beam transfers its optical vortex to 
the regenerated electric field £^^\ In the A scheme it is not allowed to have an optical 
vortex for the control beam due to adiabaticity violation at the center of the vortex. 
Using a tripod scheme for the regeneration lifts up this restriction. The probe beam 
may itself carry a vortex at the beginning of the storage [H]. Subsequently the vortex 

(s) 

is stored onto the atomic bright state $3 and then transferred back to the probe beam 
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after the control beam ^2^,2 is turned on. In that case the phase of the restored vortex 

(r) (s) (r) (s) 

in the probe beam is defined by the product f2^2 $b which both fi^2 $g may 
carry vortices. If these vortices have opposite winding numbers, they cancel each other 
leading to zero vorticity in the regenerated probe beam. 

(r) 

Let us take the restoring control laser fi^2 to be the first order Laguerre-Gaussian 
(LG) beam: fig = Ape'^ exp(— p^/o"^), where p = p/A is a dimensionless cylindrical 
radius, X = 27i/k being the optical wave-length. On the other hand, the control beam is 
assumed to be the zero-order LG beam during the storage stage involving a A system: 
exp(— p^/cTg), where a determines a relative amplitude of the control fields 

(r) (s) 

fi^2 s-iid r2^2 , '^r (Js being their dimensionless widths. This provides the following 
regenerated probe field fl?Tl) 

S^'^ = ape'^ exp[-p''{a;^ - a;^)p'^ . (32) 

3. 3. Transfer of the optical vortex during the storage of slow light 

Suppose now that initially we have a tripod system for the storage where generally both 
^^2 and are non-zero. On the other hand, a A system is used in the retrieval stage 
with only one control field, i.e. = and hence |^^2''l = 1- that case equation (!28ll 
leads to the same equation (1311) for the regenerated electric field. Yet now it is the 
storing control beam that has an optical vortex Q^^^ ~ 6^^*^. Subsequently the vortex is 
transfered to the regenerated probe field in the phase conjugated form: ~ e"'^"^. 

Suppose that the control lasers are the first and zero order LG beams at the storage 
stage: 

^2 = ApV^ exp(-pVa,^) , fig = bA eM-~pV^'s) , (33) 

where the parameter b determines the relative amplitude of the additional control laser. 
On the other hand, the control beam is assumed to be the zero-order LG beam at the 
retrieval stage involing the A scheme: fig = aAexp{—p^/a^). Thus one arrives at the 
following regenerated probe field containing the phase conjugated vortex 

£:« = ^^Pe-^^ exp[-p^(a7^ - a:')p^^ . (34) 
p^ -|- 0"^ 

3.4- Energy losses 

The proposed shemes for transferring a vortex from the control laser beam to the 
regenerated probe beam avoid non-adiabaticit (absporption) losses at the center of the 
vortex due to application of an additional control laser. Yet there are another kind of 
losses in the energy of the regenerated probe beam, because a part of the stored probe 
beam remains frozen in the medium in the form of atomic spin excitations. Let us 
estimate those losses. Suppose that the incident pulse of the probe electric field S^^^ has 
a Gaussian intensity profile with a width dp along its transversal coordinate p and a 
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Figure 2. The dependence of loses on the relative amplitude ratio between the 
control beams b calculated using equation p7p . Three circles indicate three cases 
corresponding to ap = 10 and & = 3, 10 and 30. 



length Lvac in the propagation direction z. The total electric field energy entering the 
cloud of cold atoms is 

= 27r / pdp / dz |f = -L^acl^olVJ • (35) 
Jo Jo ^ 

The electric probe field at the retrieval stage is related to the stored on in a different 
way for two cases of the storage and retrieval considered above. 

In the case of lambda storage and tripod retrieval, regenerated electric field is 
described by the equation ( l32l) . The length of the regenerated pulse is L(r). This length 
is related to the length Lvac of the initial pulse via the equation L(i.) = L^^^cVg^ /vg\ 
Since the group velocity Vg in the retrieval stage depends transversal coordinate p, the 
length L(r) is also p dependent. The total energy of electric field leaving the cloud is 

W^,) = 27r / pdp / d^l^Wp 
Jo Jo 



,2 / xexp[— 2x(crj, ^ + cr^)] 

VrLvac^O / f n / 2\ < h2 / o / 2 ' (36) 



^c3 



Here cTi-s is the width of the second control beam at the retrieval stage Q). 
6y4exp(— p^/cr^g). The energy losses may be found comparing the energies of the pulses 
before and after the interaction with the cloud. For simplicity suppose that = a^s. 
Then 



262 / 2b^\ 

exp — El -— . (37) 



p \ p / \ p 



If the intensity of the second control beam is large, b ^ ap, the ratio of the energies is 
^wMs)^cr2/(262). 
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In the case of tripod storiage and lambda retrieval the restored electric field profile 
is described by equation Nevetheless the loses are found the same as in the previous 
case ( 1371) . Several values of W(^j.)/W(^s) are given in figure O When the amplitude ratio 
of the control beams b increases, the loses are seen to increase. 

4. Concluding remarks 

We have considered propagation, storing and retrieval of slow light (probe beam) in 
a resonant atomic medium illuminated by two control laser beams of larger intensity. 
The probe and two control beams act on atoms in a tripod configuration of the light- 
matter coupling in which three hyperfine atomic ground states and one excited state 
are involved. The first control beam is allowed to have an orbital angular momentum 
(0AM). Application of the second vortex- free control laser ensures the adiabatic (lossles) 
propagation of the probe beam at the vortex core where the intensity of the first control 
laser goes to zero. Using the adiabatic approximation we have derived the equation of 
motion for the probe beam and analysed it in the case where one of the control beams 
has an optical vortex. 

Storing and release of the probe beam is accomplished by switching off and on the 
control laser beams leading to the transfer of the optical vortex from the first control 
beam to the regenerated probe field. A part of the stored probe beam remains frozen 
in the atomic cloud in a form of spin excitations, a number of which increases with 
increasing the intensity of the second control laser. We have analysed such losses in the 
regenerated probe beam and provided conditions for the optical vortex of the control 
beam to be transferred efficiently to the restored probe beam. 

We have investigated in detail two cases of storing and retrieval of optical vortices 
onto the atomic medium. In the first case the A scheme is used for the storage, whereas 
the tripod setup is exploited for the retrieval. In such a situation the vortex can be 
transferred efficiently and without a distortion from the restoring control beam to the 
regenerated probe beam. In the second case the tripod system is used for the storing 
and the A system is employed for the retrieval. The vortex is then transferred from 
the storing control beam to the regenerated probe in a phase conjugated form, so the 
regenerated probe beam acquires an opposite vorticity. The regenerated beam is then 
distorted and becomes narrower as compared to the storing beam. Thus it experiences 
a larger diffraction spreading in the subsequent propagation in the medium. 
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